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Abstract: We study parity odd transport at second order in derivative expansion 
for a non-conformal charged fluid. We see that there are 27 parity odd transport 
coefficients, of which 12 are non-vanishing in equihbrium. We use the equihbrium 
partition function method to express 7 of these in terms of the anomaly, shear viscos- 
ity, charge diffusivity and thermodynamic functions. The remaining 5 are constrained 
by 3 relations which also involve the anomaly. We derive Kubo formulae for 2 of the 
transport coefficients and show these agree with that derived from the equilibrium 
partition function. 
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1. Introduction 

Fluid dynamics is an effective description of near-equilibrium systems. Properties 
of a fluid system are always slowly varying compared to some intrinsic length scale, 
for example a mean free path. This length scale is determined by the details of 
the underlying microscopic theory. The fundamental variables in fluid dynamics are 
local velocities u^{x), temperature T{x), and all the other conserved charges or their 
chemical potentials The conservation equations for the stress tensor T'^^ and 

the other conserved currents J"^ govern the time evolution in fluid dynamics. 

The stress tensor and the conserved currents are related to the fluid variables 
{u^, T, fi"'} via constitutive relations. Since fluid systems are always slowly varying it 
is appropriate to organize the constitutive relations in terms of a derivative expansion 
of the fluid variables. At every order in the derivative expansion, the independent 
terms of the constitutive relation are constructed out of the independent derivatives 
of the fluid variables. The independent terms in the constitutive relation are mul- 
tiplied by coefficients which are functions of temperature and chemical potentials. 
These coefficients are called transport coefficients. In this paper we will study the 
transport coefficients that occur in the parity odd sector, at second order in the 
derivative expansion. These are terms constructed out of various derivatives in the 
fluid variables which are odd under parity. We will consider relativistic fluid systems 
with one additional conserved current. 

It is usually difficult to compute transport coefficients from the microscopic the- 
ory and they are generically determined from experiments. However parity odd 
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transport coefficients which occur at the first order in the derivative expansion have 
been related to quantum anomahes of the microscopic theory |jl|, H, H, ^ |^ ^. 
Our goal is to see if a similar phenomenon occurs at the second order in the derivative 
expansion of the constitutive relations. One motivation to consider second order fluid 
dynamics is that ffist order fluid dynamics is known to have problems with causality 
and numerical stability. Transport coefficients which occur at second order provide 
important constraints for spectral densities through sum rules B . Parity odd trans- 



port coefficients at second order affect the dispersion relation of chiral modes [|10 
^. This phenomenon has important experimental consequences like the spatial sepa- 
ration of particles of different chirality. In the relativistic context this phenomenon 
was ffist observed holographically |]ll| and then understood due to the presence of a 



parity odd transport coefficient at second order by |T^. 

Parity odd transport coefficients at second order has been studied earlier in [|12| 
for conformal fluids. They used the principle that parity odd terms which are even 
under time-reversal invariance should not contribute to local entropy production. 
With this principle they could constrain these transport coefficients and determine 



some of them. We will use the method developed in [|13|, |1J] to determine and 
constrain the parity odd transport coefficients. We consider non-conformal fluids 
in 3 -|- 1 dimensions which admits one anomalous charge current. This method is 
based on the requirement that the fluid equations have to be consistent with any 
equilibrium partition function. Therefore the approach flrst relies on the physical 
requirement of the existence of equilibrium. More precisely: 

• In a time independent background, that is a space-time metric with a time like 
Killing vector and background gauge flelds independent of the time direction, 
any fluid equation will admit a time independent solution. 

The second assumption is: 

• The stress tensor and the charge current evaluated on this time independent 
solution can be obtained from the partition function by varying it with respect 
to the background metric and the gauge fleld. 

This method is implemented as follows. 

1. We flrst classify all the parity odd transport coefficients till the second order 
in the derivative expansion of the stress tensor and the charge current using 
symmetries. 

2. We then evaluate the stress tensor and the charge current on the equilibrium 
fluid conflguration to the second order in the derivative expansion. 



^See g for a recent review with a complete list of references. 

^We thank Yashodhan Hatwalne for bringing this reference to our attention. 
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3. 



The equilibrium partition function is written to the second order in derivative 
expansion taking all the parity odd terms into consideration. This is also done 
based on symmetries 



4. This stress tensor and the charge current obtained from the equilibrium par- 
tition function is required to agree with that obtained from the stress tensor 
evaluated on the equilibrium fluid configuration. 



From the description of the method it is clear that only transport coefficients 
which do not vanish in the equilibrium fluid configuration will be constrained or 
determined. We will see that totally there are 27 parity odd transport coefficients. 
Out of these 12 do not vanish in the equilibrium fluid conflguration. Among the 
12, we determine 7 which we label as = 1, ■ ■ - 6 and A2. and show that they 
are related to the anomaly. The rest of the 5 are constrained by 3 relations. These 
relations are also involve the anomaly. The results are summarized in ( |1.9| ) and 
( p.. 101 ). We will then use the Kubo formula to derive two of the transport coefficients 
$1, $2 and show that it agrees with the equilibrium partition function method. The 
remaining transport coefficients seem to be related to three point functions. 

The organization of the paper is as follows. In the rest of the introduction we 
summarize our main results. In section 2 we implement the method of |jl3[ to relate 
the transport coefficients which do not vanish in equilibrium to that of the anomaly. 
In section 3 we use the Kubo formula to derive two of the transport coefficients. 



1.1 Summary of the results 

As we mentioned earlier, the aim of this note is to constrain the parity odd sec- 
ond order transport coefficients of an anomalous charged fluid in the presence of 
background electric and magnetic flelds. To deflne the transport coefficients unam- 
biguously we first must have an unambiguous deflnition of fluid variables that is the 
velocities, the temperature and the chemical potentials. {u'^,T,iJ,} We will work in 
Landau frame which is deflned by the following two conditions for the charged fluid. 

J% = -g, T^% = -^w^ (1.1) 

Let us now consider the expansion of the charge current and the stress tensor T^" 
in terms of the number of space time derivatives. This is given by 

= ^ + ^ + . . . , T'^'^ = Tj; + T^';^ + • • • , (1.2) 

where the subscript (z) refer to the number of space-time derivatives. The terms J^^.-^ 
and T^^'^ for z 7^ are all perpendicular to the velocity u^. The equations of motion 
for the fluid in presence of external electromagnetic fleld are given by the following 
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conservation laws 



V^T^" = r^-J^ + ''-je^^^'^^pR^:;,, (1.3) 
C r 

= CE,B^ + '-je-'''K^,Rl,s- 

Here = J^^^u" , B^^ = \e'^^"'^u,jj^ai3 and C is the gauge anomaly coefficient and 
Cra is the coefficient of the mixed gauge-gravitational anomaly. Note that the terms 
proportional to the mixed anomaly are fourth order in derivatives, therefore they do 
not affect the analysis of the equations of motion to 2nd order in derivatives. How- 
ever the gravitational anomaly does enter the discussion of the equilibrium partition 
function at the first order in the derivative expansion . 



We will now state the known results for the form of the stress tensor and the 
current up to first order in the derivative expansion. For i = 0, the part with 
no space-time derivative the stress tensor and current is completely determined by 
thermodynamics. At first order in derivative expansion, that is z = 1, the form of 
the current and stress tensor is explicitly known |]13[, This form is consistent 
with all physical requirements in presence of anomaly and as well as an external 



electromagnetic field |]T3|, The final result for the stress tensor and the charge 



current to first order in derivatives is given by 

[T(o)Vu ={E + P)u^u, + PG^,, 

Here E is the energy density, P the pressure and g, the charge density and G^y 
is the background metric. The variables 6, a'^'^ , V^, P and B'^ are all on-shell 
independent terms which are first order in derivatives. These are defined in table |I[ 
The variables i], (, A, and refer to the first order transport coefficients. 



Scalar s 
(1) 


Vectors 
(3) 


Pseudo Vectors 
(2) 


Tensors 
(1) 


e = v^«^ 


P^"'VyU 
V^' = {E^' - TP'^'Vui^) 


= le^"'»^s^^J^^^ 





Table 1: Data at 1st order in derivative 
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Throughout this paper, the symbol ^(^i/) on any tensor ^4^1, denotes the projected, 
traceless, symmetric part of the tensor. 

A,,., = P':p'. {^^^i±^ - ^G„,) . (1.5) 

and P'^" is the projector 

p/^'^ = u^'u'' + G^^ (1.6) 
The variable u is related to the chemical potential by 

^ = f. (1.7) 

Let us now proceed to the stress tensor and the charge current at second order 
in derivatives which we denote as [T(^2)]tiiy and J^^^y Purely from symmetry consider- 
ations, the number of independent transport coefficients upto second order is equal 
to total number of possible scalars which appear in the trace of the stress tensor, 
together with the number of possible vectors which appear in the current and the 
possible symmetric traceless tensors which appear in the traceless part of the stress 
tensor. But not all of them are independent, they can be related using the equations 
of motion given in ( |1.3|) . In Table |^ we have listed all the parity odd and on-shell 



independent scalars, vectors and tensors containing two space-time derivatives. 



Pseudo-scalars 


Pseudo Vectors 


Pseudo-tensors 


(6) 


(9) 


(12) 


Si = l^d^v 






S2 = Bf'd^u 


Vg) =e'^^"^^/.(5„J/)(9,T) 


rfj = V(^S^) 


S3 = l^d^T 






54 = B^^d^T 


Vf^) = (^i.V)S^ 


tI^J = B^i,d^)iJ 


S, = l^'V^ 






Se = Bf^V^ 


K) = 


rjS = B^^d,)T 












t'^J = B(^^Vy) 






(11) eap 



Table 2: Parity odd data at 2nd order in derivatives 
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From the table it can be seen that at second order in derivatives there are 27 
parity odd transport coefficients which appear in the current and the stress tensor. 
Therefore the most general parity odd contributions at second order in stress tensor 
and current can be parametrized as follows. 



12 



i=l 



i=l 



i=l 

Where r//j , Si and V^'^-j are defined in Table |. Our goal is to constrain the transport 
coefficients $i, Xi ^iid Aj , using the existence of the equilibrium partition function. 
It can be shown among these 27 terms only 12 can be non-zero in a time independent 
equilibrium fluid configuration. They are the first 4 in the list of scalars iSi ■ ■ ■ 1S4, the 



(1) 



ffist 2 in the list of vectors V^^^^ , V^^) the ffist 6 in the list of tensors r^y 



Therefore the analysis using the equilibrium partition function can be used to con- 
strain the 12 transport coefficients multiplying these non-vanishing terms. These are 
Xi ■ ■ ■ X4 cind Ai, A2 and $1 • • ■ $6 The final result of this analysis is the following. 



$1 = 7] hi, $2 
$5 



77 



_h dbi 



21] 62, $3 
$6 = 



= V 
2r] 



du 



$4 
db2 



2r] 



fdh\ 



T dT 



;i.9) 



A2 = 



C fdbi 2bi 



2\dT T 

_r (—- — 

X4 (, ^ 



-R2 
+ R2 



0, 



i?iTAi + 



X2 - C 



:i.io) 



0, 



where 



61 



T3 (2Cv'> 



E + P 



(E + P) 



4C9Z/ 



rp2 



21^ , ^2 



E + P 



2 1 



Note that C is the gauge anomaly coefficient and C2 is related to the coefficient of 
the mixed gauge-gravitational anomaly in ( |1.3D by 



:i.iii 
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Therefore we see that the coefficients $1, ■ ■ ■ $6 and A2 are determined in terms of 
the anomaly, the shear viscosity, the charge diffusivity A and the thermodynamic 
functions £',P, T, z/. The rest Ai,xi,---X4 are constrained by 3 relations which 
involve the anomaly. 

2. Anomalous transport from equilibrium partition function 



In sub-section we will first we briefly outline the general procedure we use to 
relate parity odd transport at the second order in derivatives to the anomaly. This 
method has been described and analyzed in detail by [|T3| and |T3|. It relies on the 



analysis of the equilibrium partition function. Thus we will be able to constrain only 
those transport coefficients which do not vanish in the equilibrium configuration. 
This section also will serve to introduce the notation and conventions used in the 



paper. In sub-section |2.2| we implement this method and derive the relations given 
in O and ([TTOl) . 

2.1 The equilibrium partition function method 

We are interested in a fluid flow on a static background metric and a static external 
electromagnetic field. The most general static metric and gauge field can be written 
in the following form. 

ds'^ = G^ydx^dx'^ = —e^^{dt + aidx^Y + gijdx^dx-' , 
Gauge Field : A/j.dx'^, 

a, ai , .Ao and A/^i are all slowly varying functions of the spatial co-ordinates {x). 
Our notations are as follows: 

• Greek indices run from 1 to 4. 

• Latin indices run from 1 to 3. 

• All Greek indices are lowered or raised by the 4 dimensional metric G^^, unless 
explicitly mentioned. 

• All Latin indices are lowered or raised by the 3 dimensional metric Qij unless 
explicitly mentioned. 

• is covariant derivative with respect to the metric G^j^i, and Vj is the covariant 
derivative with respect to the metric gij. 

• For any tensor A^^ the notation denotes the traceless symmetric part of 
the tensor, projected in direction perpendicular to the fluid velocity. 

A - papl3 f + ^Pa P'^^ A^p ^ 

^{^lv) — ^fj.^u I 2 3 

where P^i, = u^u^, + G^y is referred to as the projector. 
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Our basic assumption is that in such a background any fluid equation will admit 
a time independent solution. The stress tensor and the current, evaluated on this 
time independent solution, can be generated by varying the partition function of the 
system with respect to the background metric and gauge field. If Z is the partition 
function at temperature To then then the stress tensor and the current evaluated on 
the equilibrium are given by the following formulae. 



Toe 



2cr 



- 00 I equilibrium. 



To 



6Z 



I equilibrium 



T2 

6Z 



Stti 5Ai 



rpij 



2Tr 



equilibrium 



'Jo I equilibrium 
jil 

•J I equilibrium 



Vo 



^l^im 



bZ 



5 9 



Im 



sz_ 

T 
2f 



5Z _ ^ bZ 
6aj ^ 5Aj 



y/9 



5Z 



6g- 



Im 



(2.2) 



To 



6Z 
SAo 
6Z 



T 
\f9 



5Z 



5Z' 
5p 



5A 



where 



f = Toe-^ V 



Ao 
To 



, Ai = Ai- Aotti 



The strategy we will adopt to constrain the parity odd coefficients which occur 
at the second order is the following: 

1. Write down the most general partition function upto a given order in derivative 
expansion and consistent with all the symmetries. It will be a function of a^, 
(J and Ai and their derivatives. 

2. Vary the partition function Z, to obtain the most general possible expression 

for I equilibrium and '/'^ | equilibrium • 

3. Parametrize the most general possible fluid stress tensor and current up to some 
given order in derivative expansion using symmetries. This will give the max- 
imum number of independent transport coefficients possible constrained only 
by symmetry. We have stated the results of this analysis already in equation 
O) and Q. 

4. Evaluate the most general fluid stress tensor and the current on the equilib- 
rium solution. The final outcome will contain some of the unknown transport 
coefficients. 



5. Equate the final outcome of the previous step with the stress tensor and current 
we have already obtained by varying the partition function. 
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6. This will express the transport coefficients which appear in the stress tensor 
and the current evaluated on equilibrium, in terms of the coefficients appearing 
in the partition function. 

7. Eliminating the coefficients which appear in the partition function results in 
the desired relations among the transport coefficients. 

One might wonder that to execute the fourth step, the equilibrium solution for 
the velocity, temperature and the other fluid variables needs to be independently 
found. But as it has been explained in [13] that using this method one can pertur- 
batively determine both the solution and the transport coefficients in terms of the 
free functions appearing in the partition function. 

2.2 Partition function analysis: parity-odd sector at second order 



In this subsection we shall apply the general procedure described in to the par- 
ticular case of parity-odd sector of the charged fluid at second order in derivative 
expansion. The ffist contribution in parity-odd sector comes at ffist order in deriva- 



tive expansion. This has been analyzed in detail in section 3 of |13|. We will not 



repeat the ffist order computation here but we will extensively use their result. 
Stress tensor and current from the partition function 

Since the transport coefficients we are interested in is parity-odd sector, we will 
restrict our attention to only the parity odd part of the partition function. The 
partition function Z(2) at second order in derivatives is a gauge invariant scalar 
functional of the background metric and gauge-fields.^ Hence we need to list all 
possible parity odd scalars that can be constructed from the metric functions and 
gauge fields that contain two space-derivatives. Note that since all the functions are 
time independent no time derivatives occur. There are four such scalars: 

1. e'^'^d^ufjk, 

3. e'^'^diuFjk, 

4. e'^^difFjk, 

where z/ = ^, f = Tqc"'^, Ai = Ai-AoQi and fjk = djak-dkaj, Fjk = djAk-dkAj. 

Therefore naively, the parity odd second order partition function at two deriva- 
tives can have 4 free parameters, but two of them can be related by total derivatives. 



■^Though the system is anomalous, all the effects of anomaly i.e. the anomalous transformation 
property of the partition function under the gauge transformation is accounted for in first order 
part Z(^iy Therefore in Z(2) we need to consider only the gauge invariant scalars. 
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Ignoring the total derivative terms the most general second order partition function 
can be written as 



v(2) 



Varying this partition function with respect to the metric and the gauge fields and 
using (|2.2|) we get the following second order correction to the stress tensor and 
charge current. 



— l6 d.uF,, + To{^ 



= 2Tof 



[n(2)]*^' = 0, 



To 
T 



dT 



dMi 
'df~ 



e''%fF,, + To ) e''%ff,. 



(2.3) 



It is important to note the third equation in (|2.3| ). The fact that the spatial compo- 
nents of the equilibrium stress tensor vanish at second order will serve as an important 
constraint in determining the transport coefficients. 



Stress tensor from fluid dynamics 

We have to evaluate the fluid dynamical stress tensor on the equilibrium that is a 
time independent solution in the given static background metric and gauge field. 
This equilibrium solution for the velocity field, temperature or chemical potential in 
terms of the background metric and gauge field can also be expanded in terms of 
derivatives. We shall use the following notation. 



(2) 



T|e, = T + 5T(l) + 5T(2) + ■■ 
Z/|eq = z/ + 6u(i) + 61/(2) + ■ ■ 



(2.4) 



where Su^^.y ST(^i) and are ith corrections to the zeroth order equilibrium solution 



containing i derivatives on the background data. Now we have to substitute (2.4) in 
fiuid stress tensor and current given in (|1.4| )) and extract out the part that will be 
parity odd and which involve exactly two derivatives on the background data. This 
is the part which have to be equated with 
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From the analysis done in |jT3| we know that 



u^^ = e-^{l,0,0,0}, f = Toe-^ u = ^, 

Jo 

['^^^(i)]o = 0, 6T(^i) = 0, (5z/(i) = 0, 

[6u^i)f = ~ai[6u^^)Y, (2.5) 

[^^(1)]^ = (I) ^ + hB\ 



where 



Fjk = djAk — dkAj, 

B' = \ (e^^V,, + Ao6^^ V,.) = \^''F,, - m\ ^2 6) 

Here C is the anomaly coefficient and C2 is related to the mixed anomaly in ( [1.3|) by 

C2 = STl^Cm- (2.7) 

The presence of the mixed anomaly at ffist order in the equilibrium partition function 
was pointed out in [ P^ . 

If we also expand E, P and q in terms of a derivative expansion as 

EUg = E+5E^i) + 5Ei2) + - ■ ■ , PU = P+5P(i) + 5P(2) + - ■ ■ , q\eg = ^+^i) + ^2) + - • • 
then from (|2.5| ) it follows that 

= = (5g(i) = 0. (2.8) 

Using the fact that {u^u^ = — 1) to all order in derivative expansion we find that 

['^^(2)]o [5M(i)]*[5'U(i)]j = Parity even ~ (For our purpose). (2.9) 

Also using the Landau gauge condition on the second order stress tensor and current 
one can show that on equilibrium 

[r(2)]00 = [T(2)]^ = [J(2)]0 = 0. (2.10) 
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Using (2^), (2^), (2^) and ( 2.10 ) we get the following for the second order stress 
tensor and current evaluated on equilibrium. 



(E + P)uo[6u^2)y + 5(-2wo) + 5(CeP(5), 
6q^2)Uo + S{AVo) + 5(6/o) + S{^bBo), 



(2.11) 



where 5{—2rjafj,u), 5(C6P^,y) and (5(AV^) are the two derivative corrections of {—2r]a^u), 
{(QP^u) and (AV^) when evaluated on Su^^^y ST(^i) and Siy^i). 

Now if Q^} is a tensor which is first order in the derivative expansion satisfying 
the following two conditions 



Q^i^^u^ = at all orders and l(«'',f,p) = 0) 



(2.12) 



then one can show in general that ^ 

5g« = 5[g(i)]^ = 0. (2.13) 

Similarly if Q^f}^ is a vector which is first order in the derivative expansion satisfying 
qW^^ = at all orders and Q''^^\(u^^,f,p) = 0, (2.14) 



then it follows that ^ 



(2.15) 



sQi;' = 0. 

This argument allows us to conclude that 

di-2r]croo) = 5(CePoo) = ^^(-2^) = SiC^P^) = 6{AVo) = 0. 



Similarly S{^il^) and ^(^^5^) are the two derivative corrections of {Ci^n) and {C,bB^) 
when evaluated on 5^(1) and But and P^ are already parity odd. 

Therefore S{^il^) and 6{^bB^) are going to be parity even and hence we can set them 
to zero for our purpose. Putting all this together we have the following result for 
the components of the stress tensor and the current at the second order in derivative 
expansion 



[n(^)]oo 



^(2), 



e-6E(. 

{E + P)uo[6u^2)Y = -e''{E + P)[6u^2) 



(2.16) 









+ 5u^' 


Q, 










(1) 
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Inverting ( p.l6 ) and using ( p.3|) we can determine the second order piece of the equi- 



hbrium solution for velocity, temperature and chemical potential in terms of back- 
ground data. Now the components [H'^^)]'-? and [j*-^^]* will give rise to the constraints 
on the transport coefficients. 

We have to use ( p.l6| ) to determine SP2 and [(5^(2)]*. 

The transport coefficients are determined by demanding that ( |2.17] ) is satisfied. 



For convenience let us further split the first equation in ( ^.17[ ) in two parts, the trace 
part which is obtained by contracting the equation with gij and a traceless part which 
is obtained by subtracting the trace part of the equation from the first equation of 
( p.l7|) ^. These are given by 



Trace part : 



(2.18) 

Traceless part : 



In ( ^181) we have used to set [n(2)]^J to zero. Therefore the RHS of (|]T|) 



vanishes. 

Transport in the traceless part of the stress tensor 
In this subsection we shall analyze the second equation of ( p.l8| ). 
Traceless part : 

We now go through the following steps: 

• We first evaluate the second order contribution from a*-'. 

• We then write down the most general form for [T(2)]^^'^^ from usual symmetry 
analysis and on shell independence. In general it will contain many unknown 
transport transport coefficients to begin with. But if we restrict our attention 
to only the parity odd ones, the total number of terms is 18 ( (z = 1, ■ ■ ■ , 12) 
multiplying 12 independent traceless symmetric tensors and = 1) ■ ■ ■ 1 6) 



multiplying 6 independent scalars, see ( |1.8|) ). All these 18 terms have been 
listed in table ^ 



^Note that here trace is taken with the 3 dimensional metric. 
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One can see that all the 6 XiS are not going to appear in the combination that 
we are going to evaluate in (|2.19|) . 

So we have to evaluate the rest of the 12 terms (multiplying $jS on the equi- 
librium time-independent solution. 

For this it is sufficient to substitute only the zeroth order equilibrium solution, 
since each of these terms already contain two derivatives. 

As mentioned in section 1.2, it turns out that 6 of these twelve terms evaluate 
to zero on equilibrium. So there are only 6 transport coefficients which can be 
constrained by this equilibrium analysis. These 6 terms are the following. 

-«=V(^/.), r(2j=V(^S,>, r(3) = (a<^i.)/,), (2.20) 



r, 



fii/ 



where = t^'^'^I^UydaUp and B^ = e^^'^^UydaAp. Let us now rewrite the 
traceless part of the second order fluid stress tensor with these six terms, this 
is given by 

6 

\^{2)\^^d ~ ^alT^""^]^^ + trace part + terms that vanish on equilibrium. 

(2.21) 



a=l 



From ( [2. 191 ) we see that each of these six $a's has to be such that, they cancel 
the contribution of Sa^y when all of them are evaluated on equilibrium. 

Evaluating the spatial components of these six terms on the zeroth order equi- 
librium solution given in ( ^.5p we obtain 



r 



^(2)]ii _ _ gilg 



+ C(9^ 



{Viv)lrn + (VmZ/)/i gir 



■(VfcP/ 



(yiu)B.m + (yml^)Bi gim , „fc 



+ o(a=^ 



r 



(5)l*j 



2 3 ^ ^ 



{Vif)B^ + {Vmf)Bi gi„,.-^ 



(2.22) 



where T = -^^'^^fjk and B' = |e*-''^(Fjfc + Ao/^fc). To evaluate equation (^) 
we have extensively used section 2 of 
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The spatial components of the shear tensor gives the following contribution to 

[n(2)]ij 



- 2r] 5a'^ 



- 2ri 



[<5m(i)] 

2 3 



2r,g''g- 



il , jm 



+ 



Vi[5M(i)]„ + V„[(5M(i)]i gi^ ^ gi^ ^ 

2 ^Vfc[5M(i)] - —{dka)[6u^i)] 



(2.23) 



where 



The last term in the last line of (p.23|) results from this expansion of the pro- 
jectors. Substituting (|2.5| ) in ( p.23|) we obtain the following 



-2riSa 



-2r] 



— IT 

2 



(2)i*i 



1 fdbi 



dbi 



dT 



{5)i*i 



dbo 



-(4)lii 



db. 



T dT 



-(6)i«i 



(2.24) 



Next we have to substitute (|]23D and {^2^ in (|2]T|). Now to satisfy (|2l^) the 
coefficient of each independent expression should vanish. From examining ( p.23| ) and 
( |2.24| ) it seems that demanding every independent term to vanish results in more 
equations than the unknowns which are the transport coefficients. However from the 
structure of the equations in ( |2.19| ) it is clear that the equation admits the following 
unique solution for the transport coefficients 

'dbi 



where 



bi 



b_i dbi 



$6 = 277 



$4 
db2 



2r] 



T dT 



J3 



E + P 



2Gu^ 



-AGou 



T' 
E + P 



Gu^ 



(2.25) 



(2.26) 



and E, P, T is the energy density, pressure and temperature respectively, u = ^ refers 
to the chemical potential, and t] the shear viscosity. G is the gauge anomaly coeffi- 
cient, while G2 is related to the mixed anomaly by (|1 . 1 1| ) . The transport coefficient 
$1 gives rise to chiral dispersion relations in the shear mode ^2 . 
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Transport in the trace of the stress tensor and current 

We will now constrain the transport coefficients that occur in the trace part of the 
stress tensor and the current using the equilibrium partition function. We will first 
rewrite relevant part of the current and the trace part of the stress tensor up to 
2nd order in derivative expansion. Using (|2.3| ) and ( p. 171) we obtain the following 
equation for the trace part of the stress tensor and current at second order. 



Trace part : = - C^© + T^i-^V^f^'' + [Ti2)Y')gij 



5P2 - C5Q + 



eq 



+ 



eq 



(2.27) 



(fc) 



k=l 



T 



(2)]i 

dMi 



Current : [J^^)]' = [/'T - ^l^V] - qlSu^^)]' 
9 



dT 



e'^''{djT){dkiy)-6[AV']-q[6u^2)r 



eq 



(2.28) 



Both in (|2.27|) and (|2.28|) we have used (|1.8| ) to write the fluid stress tensor and 
the current in terms of the transport coefficients. Also in ( p. 271) to evaluate the 
combination {—2rj5a'^^ + [7'(2)]*-')fi'ii we have used (|2.24| ) and ( p.25| ). One can see that 
the transport coefficients $iS drop out. 

Using ( |2.16| ) and ( p.3p we evaluate 5i?2, ^^2 and [5w(2)]*, this results in 



5E2 = e-^'^ng) 



5q2 = e-'^jf 2 



2T^ 



dMi 



dMi 



[5u 



(2)J 



E + P 



dT 

(B'diT) - 2T 



dM2 



dT 
2T^ 
E + P 



— V- 



dM. 
~df 
dMi 

w 

dM2 



2 -dM^ 

V—r:^ 1 (/ diV), 



dT 

mT) 



— u- 



dMi 



e'^\djT){duiy). 



(2.29) 



Now we have to compute SQ and S{AV^). 



dkf 
T 



^ ( ^ - ^ ) (I'd^f) + 2b2T ) {Pdiiy) 



1 /dbi _ 2bi 

2 [df " ~f 
db2 _ &2 
df^T 



+ 



[B'd.f) + 



^ du 
dh 
du 



[&d.y). 



(2.30) 
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To derive ( p.3CI| ) we have used the following identities 



V.r = -{^) l\ = -ffd^u. (2.31) 



T ^ 

Similarly 6V^ is given by the following expression 

S{AV'} = AG'^T^M" = Ae'^''[5u^,)]^Bk = A (^|^ e'^%Bk, (2.32) 

where J^jj. = df,A,y - d^Af,. 

Now we have to evaluate the 6 scalars {Si, i = 1, - ■ ■ ,6) and the ith compo- 
nent of the 9 vectors (V^^^, k = 1, ■ ■ ■ ,9) on the equilibrium solution. We need 
the equilibrium solution only at zeroth order since all of them already contain two 
derivatives. Explicit expressions for all these terms are listed in table 

As mentioned in section 1.2, only 6 terms, 4 scalars and 2 vectors are non- 
vanishing on the zeroth order equilibrium solution. Thus the relevant parts of the 
second order current and trace of the stress tensor are given by 

[T^2)r^\trace part = P''^ [Xl{l'd,u) + X2{B^d,v) + XZ^l'd.T) + XA^B^d^T)] , 

Jf2)= Ai[e'^'^°V(aai^)(<9„T)] + A2[e^'^"Vi^„//3]. ^' ^ 

We then evaluate all these six terms on the zeroth order equilibrium solution. 

{i^d^iy) = (ra,p), 
{i^d,T) = ird,T), 

[B^^d.T) = {B%f), 



Substituting ( 2.34 ) in ( p. 33 ) we express the LHS of ( p. 33 ) in terms of the back- 



ground data and the two arbitrary functions Mi and M2 of the second order parity 
odd partition function. Now using (|2.27| ) and ( |2.3| ) we can express the 5 transport 



coeffcients appearing in ( 2.33|) in terms of Mi and M2. Since one cannot generate a 



term of the form e^^^ljBk from the partition function, the 6th transport coefficients 
A2 will be completely determined by the correction of the first order current. This is 
similar to the way the transport coefficients $t's, which appear in the traceless part 



-17- 



of the stress tensor were determined. The end result of this step is 

A hi 



A2 



2 ' 



V- 



dMi dM2 



dT dT 



X2 = -2RiT 



X3 



-2R2T 



X4 = 



dT 



'_dMi dM2 

+ c 



2\dT T 



(2.35) 



dMi\ , ^ /'db2 _ 



- fdMi 
Ai = 2T —J- I + , 

\dT \E + P 



b2 

\dT T)' 
2qT^ 



dM2 _dMi 

V 



dT 



dT 



where 



bi 



T' 



-ACou 



T' 



E + P \ 3 J ' E + P \ 2 

Ehminating Mi and M2 from these expressions we get three relations among the 
remaining 5 second order transport coefficients. 

A61 



A2 = 

T'^Ri 



T^Ri 



2 ' 



X3 - 



C fdbi 2bi 



2\dT 



T 



X4 - C 



RiTAi + 



/ db2 62 
\df " T 

- ^ Kd^ 



-R2 
+ R2 



Xi 



C(dh 
2\du 

C fdbi 



{E + P) 



2\dv 



2b2T 
-262T 



0, 



(2.36) 



0. 



3. Kubo formula for the transport coefficients $1, $2 

In this section we derive the relations obtained for the transport coefficients $1 and 
$2 given in (|2.25|) using the Kubo formula. We consider the following equilibrium 
background for the metric, gauge field and the velocity 



7(0) 



diag(-l, 1,1,1), A, 



/o^To, 0,0,0), 



;i, 0,0,0). (3.1) 



The chemical potential z/*-"^ and the temperature Tq are constants and do not depend 
on space-time. Since the energy and the pressure P^^^ can be thought of as 
functions of the temperature and the chemical potential, they are also constants in 
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space-time. Now consider the following non-zero metric perturbations about this 
background 

^Qtx = hfx, ^9tz = htzj ^9yx = hyx, Sgyz = hyz- (3-2) 
The gauge field perturbations are given by 

= (0,Oa,,0,a^). (3.3) 

The fiuid velocity is close to the rest frame and its perturbations are given by 

5u^ = {Q,v'^,Qy). (3.4) 

All perturbations are assumed to have dependence only in the time t and y-direction. 

To derive Kubo formulae for transport coefficients, we consider the constitutive 
relations for the stress tensor and the charge current as one point functions in the 
presence of external sources. We then obtain two point functions for the currents by 
differentiating with respect to the metric and the gauge field perturbations. Working 
this out to the linear order in perturbations will result in Kubo formulae for the 
transport coefficients $i,$2- 

To proceed we will require the the Christoffel symbols to the linear order in 
perturbations. The non-vanishing elements at the linear order are given by 

1 1 

Ky = -l^{dyhtx-dthyx), T\y= --{dyhu-dtKy), (3.5) 

= dthtx: ^ty — 2 i^y^tx + dthyx) , ^yy — dyh^y, 
^tt — dthzo, ^ty = 2 ^^y^^^ dthyz) , Tyy = dyh^y, 
^tx = 2 ~ dyhtx) , ^tz^ 2 ~ dyhtz) ■ 

Evaluating the inverse metric to the linear order we obtain 

h'^^htx, h'' = htz, h^^^-hyx, hy^^-hyz. (3.6) 

The covariant components of the velocity are given by 

u^^{-l,v^ + Kt,0,v' + hzt). (3.7) 

We now evaluate various components of the stress tensor to the linear order in 
the perturbations. From the list of terms that contribute at 2nd order in derivatives 
given in table 2 we see that in the background we have chosen all the scalars Si 
vanish. The reason for this is for the background all the thermodynamic functions 
are independent of space-time. We now examine the traceless part of the stress 
tensor. Note that the contributions from rjiu for i — 3,4,5,6,7,8,9,10,11 vanish 
since the thermodynamic functions are independent of space-time. What remains 
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to be evaluated are the contributions from rj^^ for i = 1,2 and i = 12. Let us first 
examine the tx and the ty component of the stress tensor. To the second order in 
derivatives and to the hnear order in the perturbation this is given by 

j^x ^ (^(0) ^ p{Oy ^ pW/i,^, (3.8) 
r*-^(E(o) + p(o))^- + p(o)/i,,, 

rpty ^ Q 

Note that a^^,a*^, [r*^*^]*^ and [r^*)]*^ for ? = 1,2, 12 do not contribute at the hnear 
order. The reason is that the term VqW/j and Vq//? is a first order term, therefore one 
has to evaluate the projector for these components say P*"P^^* at the zeroth order, 
which vanishes. 

Lets examine the yx component of the stress tensor. We need to evaluate the 
contributions from [r^*^]^^ for i — 1, 2, 12. These are given by ^ 

[r^'Y-ldliv^ + h,t), (3.9) 



We also need the contribution of the shear tensor to the linear order. This is given 

by 

a^- = i(ax + 5tM- (3-10) 

Considering all these contributions along with the contribution to the stress tensor 
to the zeroth order in derivative we obtain 

^ _p(o)/,^^ _ ^^gy + d^h^^) (3.11) 

The equations of motion for the x component of the stress tensor to the linear order 
in the fields is given by 

ajT*" + dtht^r' + dyTy^ = o. (3.12) 

Here T** is the zeroth order tt component of the stress tensor which is given by 

r« = £(°). (3.13) 

^In evaluating these contributions we take e'^^^^^ = ■y=- 
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Fourier transforming these equations and taking the zero frequency hmit or taking the 
time independent situation we obtain the Ward identities for the one point function 
of the stress tensor. 

hmT^^(A;)=0. (3.14) 

We can now differentiate these with respect to the background fields h^t, and and 
obtain Kubo formulae for the transport coefficients $i and $2 respectively. A similar 
procedure for the yz component of the stress tensor yields the same result. 

To proceed we first eliminate and using (|3.8| ). This results in the following 
equations 

r»-,-p.o.V-,( ^'^;„-fy- ) (3.15) 



2 I E(o) + P(o) / 2 
$12 f d^yT'^ - P^^^d^yht 



^(0) + P(o 



Here we have already used time-independence to drop the time derivatives. Fourier 
transforming these equations we obtain 

T-(A;) = -P(°) V - ^ ^(o) + p(o J (3-16) 
' ' — k ttz 



2 V ^(o) + P(o) / 2 



2 V + P^^^ 
= 0. 

The last equality in the above equation implements the Ward identity given in ( p.l4| ) 
and k is the momentum in the y direction. Let us now focus on the expression for 
T^^ a similar analysis can be repeated for T^^. Differentiating the Ward identity for 
T^^ with respect to hzt and and setting the other backgrounds to zero we obtain 
the following two equations 

[($i(T*^(A;)T*"(-A;)) + E^^^) + <l>i2((T*^(A:)T*"(-A;)) - p(°))] 



^(0) _^ p(0) 



<fi + $i2 ,2/^t.,,,,., ,u,^;.2_ 2r/ 



e{r%k)f{-k)) + $2^^ = -zfc— -^(T*^(A;)j-(-fc)). 



Eio) + p(o) ^ V /J V // ' ^ ^(0) _^ p(o) 

To obtain the first equation we have differentiated with respect to hf^ and set all 
the other backgrounds to zero. The second equation is obtained by differentiating 
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the Ward identity with respect to and setting the remaining backgrounds to zero. 
The equations in (|3.17|) are sufficient to determine the transport coefficients $i, $2- 



From p] we have the following results for the various two point functions. 

lim {T'%k)T'\-k)) =ik( ^(z/WrW)3 - 2C2(rW)V(°)^ , (3.18) 

/c-s>0,w-s>0 y 3 / 

lim {T'''{k)f{-k)) = ik f-(z/WT(°))2 - C2(rW)2 



These results are given in equations (79)-(81) for a system with 3 chemical potentials 
and equations (123)-(125) for a system with a single chemical potential of p. The 
definition of the variables for the two point functions used is given in equation (47). 
This reference also uses the normalization 

~^ ^ ^ C2 = ^, (3.19) 



8 327r2' 4 768^2' ^ 24' 

for the gauge anomaly and we have rewritten the chemical potential /i in terms of 
the variable ly. From |jl6| we can read out the following correlators 



lim {T'\k)r%-k)) = lim {r%k)f{-k)) = 0. (3.20) 

fc— i>0,w— i-O k-^0,u)~¥0 

These correlators are mentioned below equation (2.16) of reference [|l^. Substituting 
the formulae for the two point functions given in (|3.18|) and ( |3.20|) into the equations 
given in (|3.17|) we obtain 



*■ = + (f ("'"'r'"')' - 2C,(r<»>)».<») ) , (3.21) 

These expressions agree with that derived using the equilibrium partition function 
method which are given in ( p. 251) . 



4. Conclusions 

We have used the equilibrium partition function to obtain expressions for 7 parity 
odd transport coefficients which occur at 2nd order for a non-conformal fluid with 
a single conserved charge. These transport coefficients can be expressed in terms of 
the anomaly, shear viscosity, charge diffusivity and thermodynamic functions. Out 
of these 2 transport coefficients can also be derived using the Kubo formulae. These 
formulae agree with that obtained by the partition function method. The equilibrium 
partition function also gives 3 constraints for 5 other parity odd transport coefficients 
at this order. The transport coefficient $i affects chiral dispersion relations . 
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As we have mentioned earlier, parity odd coefficients have be examined earher 



for conformal charge transport in [|T2|. There the principle used was that these 
coefficients should not contribute to entropy production. In general the constraints 
obtained by examining the zero entropy production condition will be more than that 
obtained from the equilibrium partition function. It will be interesting to carry out 
the analysis of [jl^ to non-conformal fluids and compare with the results obtained in 
this paper. 

Our analysis of these transport properties were motived due to the possibility of 
studying them in holography. While this work was being done we received preprint 
[0 which evaluates all 2nd order transport coefficients for a charged conformal fluid 
in the framework of AdS/CFT. It is useful to compare the relations we have obtained 



for the parity odd sector with the expressions of . 

Finally it will be useful to develop Kubo like expressions for all the parity odd 
transport coefficients. From the constitutive relations it seems that the remaining 
transport coefficients involves 3 point functions. Determining these relations will 
provide an alternate method to obtain the transport coefficients from holography. 
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